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! Abstract 

■ The mass corrections to the bound state energies of the three-body muonic molecular ions 

(N; 

ppfj,,pdfi,ptfi, ddfj,, dtfj, and ftp, are determined numerically from the results of highly accurate com- 
CLf putations. The total energies and some other bound state properties of these ions are evaluated to 

very high accuracy for the bound S{L = 0) — ,P{L = 1)— and D{L = 2)— states. In these highly 

■ accurate calculations we used the most recent and accurate masses of particles mp,md,mt and 

Q ■ m^ known from high energy experiments. We also investigate some bound state properties of the 

• 1—1 ■ 

muonic molecular ions. In particular, we determine the hyperfine structure splittings of the ground 



states of the pdfi,ptp, and dtfi ions. In these calculations we used our highly accurate expectation 
values of the interparticle delta-functions obtained in recent computations. The corresponding 

in ■ 

>! hyperfine structure splittings, e.g., A12 = 1. 3400149- lO'^ MHz and A23 = 3.3518984-10'^ MHz 



for the ptfj, ion, can directly be measured in modern experiments. Analogous hyperfine structure 
splittings are evaluated to very high accuracy for all five bound S{L = 0)— states in the three 
symmetric muonic molecular ions: ppp,, ddfi and ttfi. 
PACS number(s): 36.10.+Di, 36.10.-k and 31.10.+Z. 
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I. INTRODUCTION 

In our earlier study [l| we considered the bound state spectra in the muonic molecular ions 
pp/i, pdfi, ptfi, ddfi, dtjj, and tt/i, where the notations p, d, t stand for the nuclei of three hydro- 
gen isotopes (protium, deuterium and tritium, respectively), while /i means the negatively 
charged muon u^. In our calculations in hi we have used the particle masses taken from 

n 

relatively old experimental papers, since in jl| we wanted to show the progress achieved 
recently in highly accurate computations of Coulomb three-body systems with arbitrary 
masses. Therefore, it was some logic to use the same particle masses in all such calculations 
see, e.g., 

a, a 

and references therein). On the other hand, right now the masses of all 
nuclei of hydrogen isotopes (p, d and t) and negatively charged muon /i~ are known to much 
better accuracy, than they were obtained ten years ago (in fact, eighteen years ago). It is 
clear that it is necessary to perform extensive recalculations of the bound state energies and 
other bound state properties in the six muonic molecular ions ppfi,pdfi,pt^,ddfi,dtfi and 
ttfi by using the 'recenty updated' particle masses. 

Modern highly accurate computations of the bound states in muonic molecular ions allow 
one to determine 15 - 18 correct decimal digits in the total energy E. In some cases, e.g, for 
the S{L = 0)— state in the pp/i ion such an accuracy is much higher and we can determine ~ 
21 - 22 correct decimal digits in the total energy. On the other hand, the masses of particles 
have been determined to the accuracy which corresponds to ~ 10 - 11 exact decimal digits 
only. Such uncertainties in particle masses lead to relatively large mistakes in the total 
energies and corresponding wave functions. Formally, this means an almost constant need 
of recalculation of the corresponding total energies and wave functions by solving the non- 
relativistic three-body Schrodinger equation with the new masses. Note also that in contrast 
with the two-electron atoms for three-body muonic molecular ions we cannot use various 
mass-interpolation formulas for the total energies, since they are not very accurate and 
contain not one, but two and even three different parameters (i.e. the ratios of the masses 
of particles) and none of these parameters is small. 

The main goal of this study is to perform the highly accurate computations of the bound 
states in the six muonic molecular ions ppfi,pdfi,ptfi,ddfi,dt^ and ttfi. All particle masses 
used in our calculations are taken from the most recent high energy experiments. Our 
calculations are performed with the use of extended arithmetic precision. Finally, the mass 
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corrections to the total energies of these ions have been determined (for each bound state 
in these six ions) to very high accuracy. 

By using the highly accurate expectation values of all (three) interparticle delta-functions 
obtained in our calculations we also investigate the hyperfine structure of the bound S{L = 
0)— states in the six muonic molecular ions ppfi,pdfi,ptfi,ddfi,dtfi and ttfi. The hyperfine 
structure splittings are determined for each of these (nine) bound S{L = 0)— states. The 
results of this investigation lead us to a number of interesting conclusions and observations. 
Many of these facts are important in analogous computations of more complicated systems, 
e.g., in the analysis of hyperfine structure splittings in the four-, five- and six-body quasi- 
atoms and ions which contain muonic molecular ions as a part of their structure. 



II. THE MASS-DEPENDENT HAMILTONIAN OF MUONIC MOLECULAR 
IONS 

In the non-relativistic approximation the Hamiltonian of the three-body muonic molecular 
ion abfj. (or (abfi)^) takes the form 

2mf,^ma " rUb ^' Tab r^^ 



where Vj = ^) and i = a,b, fi. In Eq.([T]) h is the reduced Planck constant {h 

and e is the elementary electric charge. It is very convenient (see below) to consider the 
bound state spectra of such ions in muon-atomic units in which h = l,m^ = 1 and e = 1. 
The speed of light c in these units is c = a~^, where a = |^ is the fine structure constant. 
In muon-atomic units the same Hamiltonian, Eq.([l]), is written in the form 

i^ = -k—V^ + -V^ + Vy + -- — -- (2) 

2^ma rUb ^' Tab Ta^j, 

where the nuclear masses rria and rrib of the two hydrogenic nuclei must be expressed in 
terms of the muon mass m^. 

In our earlier studies (see, e.g., and references therein) we used the following masses 

of the hydrogenic nuclei p, d, t and negatively charged muon fi 

= 206.768262me , mp = 1836.152701me (3) 
rrid = 3670.483014me , mt = 5496.92158me 
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where rUe designates the electron mass. In particular, these masses were used in our earlier 

n n 

studies (see, e.g., [2[, [3[ and references therein). In this work the updated values for the 
nuclear masses of all four particles p'^,d~^,t~^ and /i~ will be used. The masses of these 
four particles have recently been determined in various high-energy experiments to better 
accuracy than they were known in the middle of 1990's. Usually, these masses are expressed 
in special high-energy mass units MeV/c^. In these high-energy units the corresponding 
masses are 

rrif, = 105.65836668(38) , rup = 938.272046(21) (4) 
TUd = 1875.612859(41) , rut = 2808.290906(70) 

These values include current experimental uncertanties. They are very close to the masses 



of these particles given in 16(|. In all calculations performed in this study we have used the 
following masses of the fi,p, d and t particles (in MeV/c^) 

= 105.65836668 , rUp = 938.272046 (5) 
rud = 1875.612859 , mt = 2808.290906 

These values are considered as exact. The corresponding corrections to these masses can be 
taken into account by performing direct variational computations with the 'new' masses. 

Since our calculations of muonic molecular ions are performed in muon-atomic units, 
then we need to use only three following mass ratios — , — and — . These mass ratios 
are, in fact, the dimensionless parameters which determine the energy spectra and all other 
properties of the six muonic molecular ions ppfi, pdfi, ptfi, ddfi, dtfi and ttfi. The total energies 
and the corresponding wave functions of the muonic molecular ions are determined during 
the highly accurate solution of the non-relativistic Schrodinger equation H'${ri,r2,r^) = 

■ \['(ri, r2, rs), where E < and the non-relativistic Hamiltonian of the three-body system 
which is written in the form of Eq.(|2]). To determine the highly accurate solutions of the 
non-relativistic Schrodinger equation with E < in this study we apply the exponential 
variational expansion in the relative/perimetric three-body coordinates. The explicit form 
of the exponential variational expansion in perimetric coordinates is 

N 



'^LM = ^(1 + /«Al) '^i^LM'(r31, r32) exp{-aiUi - PiU2 - ^iU^) (6) 



'-1 
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where Ci are the hnear (or variational) parameters, ai,(3i,'yi are the non-linear parameters 
and L is the angular momentum of the three-body system abn. Note that each basis function 
in Eq.([H]) is an eigenfunction of the and operators with eigenvalues L{L + 1) and M. 
This means that LF'^! = L{L + while M is the eigenvalue of the operator, 

i.e. Lz'^ LM = M^! LM- The operator P21 in Eq-® is the permutation of the two identical 
particles in symmetric three-body systems. For such systems in Eq.([6]) one finds n = ±1, 
otherwise k = 0. In general, for the bound states of natural spatial parity we chose in 
Eq.Q K, = (—1)^ for all symmetric muonic molecular ions ppfi,ddfi,tt^ and k = for all 
non-symmetric ions pdfi,ptfi, dtfi. ^ 

The functions 3^^^j{^^(r3i, r32) in Eq. fll4p are the bipolar harmonics [5] of the two vectors 
Tsi = '"si ■ n3i and r32 = r32 ■ n32. The bipolar harmonics are defined as follows sl 

where C^^^.^^^^ are the Clebsch-Gordan coefficients (see, e.g., [5| and [6|) and the vectors 
Ha; = f and iiy = ^ are the corresponding unit vectors constructed for arbitrary non-zero 
vectors x and y. As follows from Eq.(I7]) each bipolar harmonic is the M— component of 
the irreducible tensor of rank L. In actual calculations it is possible to use only those 
bipolar harmonics for which £1 + £2 = L. Note that the basis set, Eq.(l6]), is a partial case 
of the more general exponential variational expansion in the relative/perimetric coordinates 
In particular, our Eq. (fT4|) does not include exponents with the imaginary (or complex) 
non-linear parameters and some other factors which are needed to accelerate the overall 
convergence rate for some three-body systems, e.g., for the Ht, Dj" ions, helium-muonic 
atoms and for other 'special' systems (for more details, see, 1| and j?!). We also do not 
want to discuss here the bound states of unnatural spatial parity, when one needs to use in 
Eq.dnj the bipolar harmonics for which ^l + ^2 = L + 1. 



III. SPECTRA OF BOUND STATES IN MUONIC MOLECULAR IONS 

By analyzing the bound state spectra in the six muonic molecular ions 
ppfi,pdu^tfi,ddfi,dtfi and tt^ one finds that they can be separated into three different 
groups 2| on qualitative grounds. The first group includes three light muonic molecular 
ions ppfj,,pdfj, and ptfi. Each of these systems has two bound states: one S{L = 0)— state 
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and one P{L = 1)— state. Neither of these two states is weakly bound. Note that each of 
these hght muonic molecular ions contains at least one protium nucleus. The second group 
includes the two 'intermediate' muonic molecular ions ddfi and dtfi each of which has five 
bound states: two S{L = 0)— states, two P{L = 1)— states and one D{L = 2)— state. One 
of these five states (the excited P*{L = 1)— state in each of these two ions) is weakly bound. 
By the formal definition the weakly bound state in a few-body system is a state with very 
small binding energy e, or, in other words, with a very small ratio of the binding and total 
energies, i.e., r = -| ^ 1, where r is the dimensionless parameter. The third group con- 
tains only the heaviest muonic molecular ion tt/i which has six bound states (and no weakly 
bound states): two S{L = 0)— states, two P{L = 1)— states, one D{L = 2)— state and one 
F{L = 3) -state. 

The classification of bound state spectra in muonic molecular ions is based on the general 
theory developed in j^, {q] for three-body Coulomb systems with unit charges. This theory is 
based on the fact that the total number of bound states in any muonic molecular ion a^b^^~ 
is determined by the lightest nucleus of the hydrogen isotope in this ion. This explains why 
only three groups of different bound state spectra can be found in the six muonic molecular 
ions mentioned above: the p— group, the rf— group and the t— group. Furthermore, it must 
be a similarity between the spectra of bound states in each group: e.g., between the bound 
state spectra of the 'protium' muonic molecular ions ppfi,pdfi and ptfi. Analogous similarity 
can be found for the bound state spectra of the ddfi and dtfi ions in the 'deuterium' group. 
It can be shown that in such 'families' of muonic molecular ions the symmetric ion always 
has the maximal binding energy jo]. By using these similarities between the bound state 
spectra in each of these 'families', one also finds a number of useful relations for the total and 
binding energies as well as for other bound state properties of different muonic molecular 
ions (see examples in [l|). 

As we have mentioned above there are 22 bound states in the six muonic molecular ions 
ppfi,pdfi,ptfi, ddfi, dtfi and ttfi. Nine of these states are the S{L = 0)— states, while nine 
others are the P{L = 1)— states. There are also three bound D{L = 2)— states and one 
bound F{L = 3)— state. The F— state is stable only in the heavy ttfi ion. In this study we 
determine the total energies of twenty one such bound states. At this moment we cannot 
perform the highly accurate computations of the F{L = 3)— state in the ttfi ion, since our 
unique code for such calculations was lost a few months ago (due to some problems at our 
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local computer). 



IV. RESULTS AND DISCUSSIONS 

In this study all numerical computations of the S— and P— bound states in muonic 
molecular ions in this study are performed with the use of 64 - 108 decimal digits per 
computer word 10|, ll|], allowing the total energies to be determined to the accuracy 
^ 1 • 10^^^ — 1 • 10~^^ m.a.u. In all calculations we have assumed that all particle masses 
and corresponding conversion factors (e.g., the factor Ry mentioned below) are exact. In 



fact, such assumptions are always made in the papers on highly accurate computations 



in few-body systems (see, e.g., 12| and 13|). The known experimental uncertainties in 



particle masses and conversion factors are taken into account at the last step of calculations, 
when the most accurate computations are simply repeated for a few times with the use 
of different particle masses and conversion factors. Analogously, the expectation values of 
other operators are determined in calculations with our non-relativistic wave functions. To 
avoid a substantial loss of numerical accuracy during computations of the expectation values 
of some operators these non-relativistic wave functions must be extremely accurate. 

Table I contains the total variational energies obtained for the ground and first 'vibra- 
tionally' excited S{L = 0)— states of the symmetric and non-symmetric muonic molecular 
ions pp^,pdfi,ptfi, ddfi, dtjj, and tt/i. Table II includes the total energies for the rotationally 
and vibrarionally excited P{L = 1)— and P*{L = 1)— states of these six muonic molecular 
ions. In these two Tables and everywhere below the upper index is used to designate the 
vibrationally excited state with the same angular momentum L. Table III contains the total 
energies of the D{L = 2)— states in the three heavy muonic molecular ion ddfi,ttfi and dtfi. 
The results from Table III have been obtained with the use of the standard Fortran with 
the quadruple precision accuracy (30 decimal digits per computer word). Also, in calcula- 
tions of the bound D{L = 2)— states we did not applied our two-stage optimization of the 
non-linear parameters in the trial wave functions. Therefore, the energies from Table III are 
less accurate than analogous energies from Tables I and II. 

In Tables I - III the notation oo stands for the total energy which corresponds to the 
infinite number of basis function, i.e. = oo in Eg. (1141) . The asymptotic formula for the 
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total energy takes the following four-parameteric form 

E{N,) = Eioo) + ^ + (8) 

where E{oo), Ai, A2 and 7 are the four parameters which are determined by using the results 
of highly accurate calculations of the total energies with the different numbers of basis 
functions Ni (see Eq.([6])). To determine four parameters in Eq.® one needs to know at 
least four total energies E{Ni) obtained from the direct numerical calculations. In general, 
the asymptotic value of the total energy contains one/two correct decimal digit (s) extra. 
Moreover, the exact coincidence of some decimal digits in the E{Ni) and E{oo) energies in 
Eq.(|8]) allow us to confirm the total number of stable decimal digits in the final energies. 

Highly accurate calculations of the total energies and other bound state properties for all 
known bound states in muonic molecular ions is an important and actual scientific problem. 
As mentioned above in this study we use the improved values of particle masses known from 
recent high energy experiments. Our results presented in Tables I - III provide answers for 
a number of actual questions. For instance, as is well known the bound P*{L = 1)— states 
in the ddfi and dtfi ions are very weakly bound. By using the corresponding energies from 
Table II and improved masses of the muon, deuterium and tritium nuclei (see above) one 
finds the 'improved' binding energies of the P*{L = 1)— states in the ddfi and dtfi ions: 

£{ddfi; P*{L = 1)) = -1.9749828376301(2)el^ (9) 

and 

e{dtfx; P*{L = 1)) = -0.6603325645(2)e1/ (10) 

where we have used the following conversion factor Ry = 27.211385060 ■ (^) from muon- 
atomic units to electron vohs (1 eV = 1.602176487(40) xlO"^^ J). Note that such a re- 
calculation from muon-atomic units to electron volts requires the knowledge of the electron 
mass rrie (me = 0.510 998 910 MeV/c^) and Rydberg constant Ry (or conversion factor). 
The binding energies of any other bound state in muonic molecular ions can be evaluated 
analogously. These values can be compared with the binding energies e{ddfi; P*{L = 1)) 
and €{dtix; P*{L = 1)) determined in [l|. 

Another interesting problem is to study the changes in the bound state properties of these 
muonic molecular ions which are directly related to the mass variations. The results of these 
calculations (in muon-atomic units) can be found in Table IV for some of the properties. 
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The results from this Table can be compared with analogous results from Table 8 given in 
Such a comparison shows the effect of mass variation for the bound states properties 
which are different from the total and binding energies. By working with Table IV we have 
found the numerical mistake in the (r|^) expectation value for the pdfi ion (extra '1' was 
added from the left side to this result in Table 8 in The mass corrections to the total 
energies can be obtained by subtracting our results from Table I - III from the corresponding 



results given in Tables 1 - 5 in 



l|. Since the total energies in all these Tables are given in 



muon-atomic units, then we do not need to use any additional conversion factor. 

V. HYPERFINE STRUCTURE OF THE GROUND STATES IN THE pdfi,ptn AND 
dtn IONS 

In this Section we analyze the hyperfine structure and determine the hyperfine structure 
splitting of the bound S{L = 0)-states in the non-symmetric muonic molecular ions pdfi,ptfj; 
and dtfi. As is mentioned above there are four bound S{L = 0)— states in these three ions: 
three ground S{L = 0)— states (one in each of these ions) and one excited S{L = 0)— state 
in the heavy dtfi ion. In this Section we want to investigate the hyperfine structure and 
determine the hyperfine structure splitting for each of these bound states by using highly 
accurate expectation values of the delta-functions obtained in our highly accurate numerical 
computations (see above). 

The general formula for the hyperfine structure splitting {AH)h s (or hyperfine splitting, 
for short) in the case of an arbitrary three-body system is written as the sum of the three 
following terms. Each of these terms is proportional to the product of the factor and 
expectation value of the corresponding (interparticle) delta-funtion. The third (additional) 
factor contains the corresponding (7— factors (or hyromagnetic ratios) and scalar product of 
the two spin vectors. For instance, for the pdfi ion this formula takes the form (in atomic 
units) (see, e.g., Cl) 



(Ai/)... = ^a2^(5(r,,))(s, ■ s,) + ^a^^^(5(r,,))(s, ■ s,) 

(11) 

2 

where a = |^ is the fine structure constant, and nip are the muon and proton masses, 
respectively. The factors g^^Qp and Qd are the corresponding (7— factors. The expression 



for {AH)h.s. is, in fact, an operator in the total spin space which has the dimension (2sp + 
l)(2sd + l)(2s^ + 1) = 12. In our calculations we have used the following numerical values for 
the constants and factors from Eq.f llip : a = 7.297352586 ■ 10~^, nip = 1836.152701me, = 
206.768262me and = —2.0023218396. The (7— factors for the proton and deuteron are 
deteremined from the formulas: Qp = ^ and = where Aip = 2.792847386 and 
Aid = 0.857438230 are the magnetic moments (in nuclear magnetons) of the proton and 
deuteron, respectively. The spin of the proton and deuteron is designated in Eq. ffTT]) as 
/p = i and Id = 1. 

The analogous formula for the hyperfine structure splitting in the ptn ion takes the form 

3 3 mpTn^ 

where gt = where Ait = 2.9789624775 is the magnetic moment of the triton expressed 
in the nuclear magnetons and It = ^ is the spin of the triton (or tritium nucleus). The 
formula for the hyperfine structure splitting in the dtfi ion is 

(AH),.,. = ^a'^{6irdt))isd ■ s,) + ^a2^^(5(r,,))(s, • s,) 
3 nip 3 mpm^ 

+ ^«^^(5(r*,))(s,-s,) (13) 
3 mpTrin 

where all values are defined above. The same formula can be applied to determine the 
hyperfine structure spllitting in the excited S{L = 0)— state of the dt^ ion. The only 
difference in the hyperfine structure splittings determined for the ground and excited states 
of the dtfi ion can be related with the expectation values of interparticle delta-functions. 

In our computations of the muonic molecular ions performed recently |l| we have deter- 
mined the expectation values of all delta-functions which are needed in Eqs.f llip - f|T3|) . The 
corresponding expectation values are shown in Table I. These values have been determined 
in muon atomic units where = l,h = l,e = 1. They must be re-calculated to the regular 
atomic units (mg = l,h = l,e = 1) which are used in the formulas, Eqs. fllip - f ll3p . to 
determine the hyperfine structure splittings. In these calculations we have used the trial 
wave functions with N = 3300, 3500, 3700 and 3840 exponential basis functions (for more 

n 

details, see [l[). The expectation values of all interparticle delta- functions computed for the 
ground S{L = 0)— state of the pd/j, ion are shown in Table I. The overall convergence rates of 
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the delta-functions computed for each bound state in the ptjj and dtjj, ions are very similar 
to the results shown in Table V. 

These expectation values of the S{rij) functions were used in the formulas Eqs. fllip - 
f[T^ to determine the hyperfine structure splittings of the bound S{L = 0)— states of the 
pd^,ptfi and dtfi ions. Numerical values of the corresponding hyperfine structure splittings 
can be found in Tables VI and VII. Note that these values are usually given in MHz, 
while the values of {AH)h,s. which follow from Eqs. ffTTl) - (1131) are expressed in atomic units. 
To re-calculate them from atomic units to MHz the conversion factor 6.57968392061 -10^ 
MHz/a.u. was used |l6 |. 

In general, the pdfi and dtfi ions have similar hyperfine structure. In particular, in each 
of these ions one finds twelve spin states which are separated in the four following groups: 
(1) the group with J = 2 (five states), (2) the group with J = 1 (three states), (3) the 
group of one state with J = (one state) and (4) the group with J = 1 (three states). 
Here and everywhere below the notation J stands for the total spin (or total momentum, 
for the S{L = 0)— states) of the three-body ion. The states with J = 2 have the maximal 
energy, while the energy of the states from the fourth gourp is minimal. The corresponding 
splittings Ai2, A23 and A34 can be found in Table VI for each bound state in the pdfi and 
dtfi ions. 

The hyperfine structure of the ground state in the ptfi ion is completely different (see 
Table VII), since the spin of the triton equals |, while the spin of the deuteron (or deuterium 
nucleus) equals 1. In the case of the ground state in the ptfi ion one finds only eight spin 
states which are separated into three different groups: (1) the group of four states with 
J = |, (2) the group of two states with J = \ and (3) the group of two states with J = |. 
The group (1) has the maximal energy, while the energy of the states from the third group is 
minimal. The corresponding values of the hyperfine structure splittings in the ground state 
of the ptfi ion are A12 = 1.3400149-10^ MHz and A23 = 3.3518984-10^ MHz. 

In this Section we have investigated the hyperfine structure and determine the hyperfine 
structure splitting in the bound S{L = 0)— states of the pdfi,ptfi and dtfi ions. The first 
excited S{L = 0)— state in the dtfi ion is traditionally designated by an additional asterisk, 
i.e. (dtfi)*. In such calculations we used the highly accurate expectation values of all inter- 
particle delta-functions obtained in recent computations [l|. In general, it is very interesting 
to compare the numerical values of the hyperfine structure splittings A12, A23 and A34 for 
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different muonic ions (see Table VI). 



VI. HYPERFINE STRUCTURE OF THE BOUND S{L = 0)-STATES IN THE 
SYMMETRIC MUONIC MOLECULAR IONS 

In this Section we consider the hyperfine structure splitting in the symmetric muonic 
molecular ions ppfi, dd^ and tt/i. As is well known there are five bound S{L = 0)— states in 
these (symmetric) muonic molecular ions. The ground states are stable in each of these ions, 
while the excited S{L = 0)— states are stable only in the heavy ddfi and ttfi ions. In general, 
the analysis of the hyperfine structure in symmetric systems is slightly more complicated 
than analogous analysis for non-symmetric systems/ions. On the other hand, the arising 
hyperfine structure is relatively simple and can be explained by using a few transparent 
physical ideas. 

The general formula for the hyperfine structure splitting (or hyperfine splitting, for short) 
for an arbitrary three-body muonic molecular ion aa/i is written in the following form (in 
atomic units) (see, e.g., [14]) 

(AH),,,, = ^«2^(5(r,J)(s. ■ sj + ^a'^^{S{r,,)){s, ■ s,) 

+ ^«^^^(5(M)(s„-s,) (14) 

2 

where a = |^ is the fine structure constant, and rup are the muon and proton masses, 
respectively. The factors g^^ and Qa are the corresponding factors. The expression for 
{AH)hs^ is, in fact, an operator in the total spin space which has the dimension {2s a + 
1)2 (2s^ + 1). Since the second and third terms in Eq.(l) are identical, then we can reduce 
Eq.dUl) to the form 

{AH),.,. = '^a'^-^{6{r^a)){sa ■ s^) + '^a'-^{6{r^,)){Saa ■ s,) (15) 
o Trip 6 mpm^ 

where Saa = (sa + Sa) is the total spin of the pair of identical particles (the two nuclei of the 
hydrogen isotopes), i.e. p,d and t. 

The formula, Eq. ffTSl) . allows one to make a few qualitative predictions about the hyperfine 
structure of the symmetric muonic molecular ions. First, it is clear that the classifications 
of the levels of hyperfine structure must be based on the total spin of the two 'symmetric' 
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nuclei Saa- The absolute values of the spin Saa are always non-negative integer numbers, 
i.e. I Saa 1= 0, 1, 2, . . .. For instance, in the case of two protons p and/or two tritons t one 
finds I Saa 1= 0, 1, while for the two deuterons | Saa |= 0, 1, 2. The hyperfine energy of this 
state with J = is determined only by the first term in Eq. flT^ which is very small, since 
the expectation values (5(raa)) in all muonic molecular ions are very small. As follows from 
actual computations all these values are less than 4- 10~^ (in muon atomic units). Briefiy, we 
can say that the energy of this hyperfine state (with J = i) is determined by the spin-spin 
interaction between the two heavy nuclei (muon's spin does not contribute). The overall 
contribution from the first term in Eq. f[T^ rapidly (exponentially) decreases when the mass 
of the heavy particle increases. Formally, the first term in Eq. fllSp is very small already for 
the ppfi ion. However, for the ddfi and ttfi ions its contibution is negligible. This means 
that in the first approximation the hyperfine structure of the symmetric muonic molecular 
ions can be explained by using only one term for the muon-nuclear spin interaction. This 
leads to some 'additional' symmetry observed for the actual levels of hyperfine structure of 
heavy ions (see below). 

As is well known the spin of the negatively charged muon /i~ equals | and the spins of 
the proton p and triton t also equal |. Therefore, the hyperfine structure of the ppfi and 
tt/i ions must include eight levels which form three following groups: (1) the group of four 
spin states with J = |, (2) the upper group of two states with </ = | and (3) the lower 
group of two states with J = |. The classification is true for the excited S{L = 0)— state in 
the ttfi ion. Here and everywhere below the notation J stands for the total spin (or total 
momentum) of the three-body system/ion, since J = L + S = S for the S{L = 0)— states). 

The hyperfine structure of the ddfi ion is substantially different. In the ddfi ion one finds 
eighteen levels of hyperfine structure which are separated into five different groups: one 
group with ^ = I (six states), two different groups of states (upper and lower groups) with 
J = I (four states in each), two different groups of states (upper and lower groups) with 
J = i (two states in each). 

In our calculations of the hyperfine structure we have used the following numerical values 
for the constants and factors in Eq.([l5]): a = 7.297352586 ■ lO^^,^^ = -2.0023218396 
and rUp = 1836.152701me, = 206.768262me. The (7— factors for the proton, deuteron 
and triton are deteremined from the formulas: Qp = ^j^,gd = and gt = where 
Mp = 2.792847386, Md = 0.857438230 and Mt = 2.97896247745 are the magnetic moments 
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(in nuclear magnetons) of the proton, deuteron and triton, respectively. Here the spins of 
the proton, deuteron and triton are designated in by the letter I with the corresponding 
index: Ip = ^, 1^ = I and = |. In Eqs. lini) - ( IT^ these values are designated differently. 
In highly accurate computations of the expectation values of delta-functions we have used 
the following masses of the deuteron and triton: = 3680.483014 rrie and rrit = 5496.92158 
rrie. These masses are often used in modern highly accurate calculations of muonic molecular 
ions (see, e.g., [l|). 

The convergence of the expectation values of the nuclear- nuclear (or pp—) and nuclear- 
muonic (or p/i— ) delta-functions is illustrated in Table VIII for the ppfi ion. The convergence 
of these expectation values for other bound S{L = 0)— states in the ddfi and ttfi ions is very 
similar to the results presented in Table VIII for the ppfi ion. The hyperfine structure and 
energy splittings between the corresponding levels for all five bound S{L = 0)— states in the 
three muonic molecular ions ppfi, ddfi and ttfi can be found in Tables IX and X. In atomic 
physics these values are traditionally given in MHz. The corresponding conversion factor is 
6.57968392061-10^ MHz/a.u. In Tables IX and X the excited states are designated by the 
asterisk used as the upper index, e.g., (ddfi)* and (ttfi)*. Such a system of notation is often 
used for muonic molecular ions. 

Tables IX and X contain both the energies of the levels of hyperfine structure (ej) and 
hyperfine structure splitting (A(Ji — t- J2)). As we have predicted (see above) one of the 
hyperfine levels has a very small energy. As follows from Tables IX and X this level cor- 
responds to J = |. In the ttfi ion the hyperfine energies of this state are ~ 11.0591 
MHz and ~ 12.4307 MHz for the ground and first excited states, respectively. In the 
ddfi ion the energies of the analogous levels are 6.8996 MHz and 4.7378 MHz, respectively. 
Briefly, this means that the overall contribution of the nuclear-nuclear spin interaction is 
very small for the both ddfi and ttfi ions. This directly follows from the known fact (see, e.g.. 



171]) that the expectation values of nuclear- nuclear delta-functions are very small. For in- 



stance, for the ddfi and ttfi ions the expectation values of nuclear-nuclear delta-functions are 
{Sdd) ~ 2.43871205-10-6 (m.a.u.), {Sdd) ~ 1.67460229-10"^ (m.a.u.), {Su) ~ 2.15893994-10^^ 
ijn.a.u.) and {5tt) ~ 2.42670033 - 10"'' (m.a.u.), for the ground and excited states, respec- 
tively. Finally, the observed hyperfine structure of these two ions is mainly (99.9999 %) 
related to the muon-nuclear spin interactions only. In the ppfi ion the situation is slightly 
different, but even for this ion the overall contribution of the muon-nuclear spin interaction(s) 
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is substantially larger than the contribution from the nuclear-nuclear spin interaction. 



VII. CONCLUSION 

Thus, we have determined the total energies of the twenty one bound states in the six 
muonic molecular ions ppfi,pdfi,ptfi,ddfi,dtfi and ttfi. The angular momentum of these 
bound states equals L = 0, L = 1 and L = 2. Our calculations of the bound state energies 
in this study have been performed with the improved particle masses known from recent 
high-energy experiments. The highly accurate wave functions of muonic molecular ions are 
needed to determine the expectation values of some operators. These expectation values 
are of interest in numerous applications related to the muonic molecular ions. Currently, all 
muonic molecular ions can be created in real experiments and their various properties can 
be measured to very good accuracy. Therefore, we can compare the predicted (or computed) 
values of the bound state properties with their actual (or observed) values. In general, the 
analysis of these three-body systems is significantly more interesting and informative than 
the traditional analysis of the two-electron atoms and ions. 

The hyperfine structure of all nine bound S{L = 0)— states of the ppfi,pdfi,ptfi, ddfi, dtfi 
and ttfi has also been investigated. In our calculations of hyperfine structure and hyperfine 
structure splittings for each of these states we used the highly accurate expectation values 
of the interparticle delta-functions. The hyperfine structure splittings of the ground states 
of each of these ions (see Tables VI, VII, IX and X) can directly be measured in modern 
experiements. 

Appendix 

In this Appendix we briefly discuss the history of the bound state computations of muonic 
molecular ions. Note that the interest to these ions was always closely related to the problems 
of muon- catalyzed fusion of the nuclear reactions. The first numerical computations of the 
30und states in three-bodymuonic molecular ions were performed by Belyaev et al in 1959 
isl . By using a very simple adiabatic (but non-variational!) procedure they were able to find 
20 bound states in six ions ppfi, pdfi, ptfi, ddfi, dtfi and ttfi. Unfortunately, due to lack of good 



computers at that time the overall accuracy of the procedure used in 18| was very low and the 



authors could not confirm the boundness of the excited P*{L = 1)— states (or (l,l)-states) in 
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the ddn and dt/j, ions. It was concluded only that, if such states are bound, then they are very 
weakly bound. The binding energy of these two states was expected to be smaller than 4.5 
eV, i.e. smaller than the binding energy of a typical hydrogen molecule. Immediately after 
publication of [isl an intense stream of speculations started about a possible interference 
(or resonance) between the formation of excited P*{L = 1)— states (or (l,l)-states) in the 
ddfi and dt^ muonic molecular ions and different transitions in surrounding molecules (see, 
e.g., [19I and references there in). 

In the middle of 1960's Halpern [20|, Carter [2l|, [22| and Delves and Kalotas [23| pub- 
lished their papers with the results of variational computations obtained for some bound 
states in the muonic molecular ions. In particular, Halpern 20| considered the bound 
P{L = 1)— states in the symmetric ppfi,ddfi,ttfi ions. The paper of Delves and Kalotas 
has a great methodological interest, since 99 % of all modern methods for highly accu- 
rate computations of Coulomb three-body systems are based on that work. At the same 
time a number of experiments have been performed by Bystritskii et al (see 24 1 and {25 1 
and references therein). They worked with /x~-muons which were slowing down in liquid 
deuterium and deuterium-tritium mixture. Finally, it was found that one muon can cat- 
alyze approximately 10 - 20 (d, c/)— nuclear reactions in liquid deuterium (D2) and 90 - 110 
((i, t)— reactions in the liquid equimolar deuterium-tritium mixture (D2 : T2 = 1:1). Such 
very large numbers of nuclear reactions catalyzed by one muon can be explained only by the 
resonance (or very fast) formation of ddfi and dtfi muonic molecular ions. Correspondingly, 
the related processes were called 'resonance' muon- catalyzed fusion of nuclear reactions, in 
contrast with the 'regular' muon-catalyzed fusion observed in 26 1. 

Those experimental works produced a great interest to study the weakly-bound states 
in the ddfi and dtfi ions. The main goal of all following computations was to determine 
the binding energies of the weakly-bound P*{L = 1)— states in the ddfi and dtfi ions to the 
accuracy ^ 1 ■ 10~^ eV (or approximately 10 K). At that time a large number of bound state 



computations for muonic mo 



ecular ions were performed with the use of the adiabatic (but 

bound 5YL = 0)- 



non-variational!) method |27|. The first variational computations of al 
and P{L = 1)— states in muonic molecular ions have been conducted in 28| and 29|. Later, 
we have substantially improved the accuracy of such computations jsO] and were able to 
calculate the bound D{L = 2)— states in these ions [sij]. Since then we have performed 
a number of different computations of the bound D{L = 2)— states in the ddfi,dtfi and 
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tt/i ion. The total energy of the dt/j, ion [32] is the only result obtained for these states 
outside of our group. However, in the middle of 1980's the masses of all particles involved 
in the muonic molecular ions were determined to much better accuracy. In addition to this, 
it became finally clear that the main restriction of the resonance muon-catalized fusion is 
directly related to the muon stiking coefficient to the ^He nucleus, or with its inverse value 
which equals to the number of nuclear reactions catalyzed by one muon in the equimolar 
deuterium-tritium mixture. It appears that such a number (150 - 160) was very close to the 



value obtained earlier 



2J] and 25| . It was also shown that even 200 nuclear fusions per one 



muon in the equimolar deuterium-tritium (liquid) mixture is in 12 - 15 times less than it 
is needed for theoretical break-even and ~ 55 - 65 times smaller than necessary for actual 



break-even (see discussion and references in 



33|). 



After these publications the overall interest to the resonance muon- catalyzed fusion 
rapidly went down. Nevertheless, in the middle of 1990's we have performed a series of 
highly accurate computations of the muonic molecular ions 3^, 351]. These works were 
originally stimulated by the development of the theory of bound states in the Coulomb 
three-body systems with unit charges. Later, I have wrote the code to determine the en- 
ergy and some other properties of the F{L = 3)— state in the ttfi ion. In our calculations 
performed around 2001 - 2003 we have used the advanced Fortran pre-translator written by 
D.H. BaUey Q, Q. Abound that t.n,e anothe. pape. was published on highly accurate 
bound state computations of the three-body muonic molecular ions [38[ . In [l[ we have used 
very large basis sets and performed an accurate optimization of the non-llinear parameters 
of our method. The paper |l| contains the most accurate values of the total (and binding 
energies) of all 22 bound states in the six muonic molecular ions pp^,pdfi,ptfi, ddfi, dtjj, and 
tt^i. In this study we wanted to recalculate some of these systems by using the new values 
of the particle masses. 

Note in conclusion that despite an obvious failure of the resonance muon-catalyzed fusion 
in the D:T mixture the muon-catalysis of the nuclear reactions is not a closed problem. In 
this area one still finds dozens of unsolved, approximately and wrongly solved problems. One 
of such problems is the current (very large) deviation between theoretical and experimental 
fusion rates in the pt/i ion. Other open problems include accurate evaluations of the fusion 
rates and muon stiking probabilities for some rotationally excited states, probabilities of 
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excitations (and de-excitations) of muon-molecular ions, etc. 
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TABLE I: The total energies (E) of the bound S{L = 0)— states in the symmetric muonic molecular 
ions in muon-atomic units (m^ = l,^ = l,e = l). N designates the number of basis functions 
used in Eq.(2). 



N 




E{ddn) 


E{ttn) 


3300 


-0.494386 815212 835026 521839 


-0.531111 132193 187917 45 


-0.546374 225613 816728 844 


3500 


-0.494386 815212 835026 522038 


-0.531111 135402 386449 61 


-0.546374 225613 816728 849 


3700 


-0.494386 815212 835026 522184 


-0.531111 135402 386450 59 


-0.546374 225613 816728 855 


3840 


-0.494386 815212 835026 522266 


-0.531111 135402 386451 22 


-0.546374 225613 816728 856 


oo 


-0.494386 815212 835026 52250(10) 


-0.531111 135402 386455(2) 


-0.546374 225613 816728 90(3) 


N 


E{pdiJ,) 


E{ptn) 


E{dtii) 


3300 


-0.512 711 792 481 703 484 


-0.519 880 085 704 058 459 


-0.538 594 971 709 480 710 


3500 


-0.512 711 792 481 703 573 


-0.519 880 085 704 058 570 


-0.538 594 971 709 480 718 


3700 


-0.512 711 792 481 703 647 


-0.519 880 085 704 058 670 


-0.538 594 971 709 480 724 


OO lU 


-U.olz (11 (92 181 (().>()(() 


-(J.-)19 880 U8;) (Ul U-)8 (11 


n -OO -n 1 n"^ 1 "^nn ion Ton 

-U.-).>8 091 9(1 (U9 180 (.j() 


OO 


-0.512 711 792 481 703 85(4) 


-0.519 880 085 704 058 87(4) 


-0.538 594 971 709 480 79(2) 


N 


E*{dd^i) 


E*{dt^i) 


E*{ttii) 


3300 


-0.47970 63771 01901 40596 


-0.488 065 354 215 765 737 


-0.49676 28898 97946 30625 


3500 


-0.47970 63771 01901 40638 


-0.488 065 354 215 765 800 


-0.49676 28898 97946 30709 


3700 


-0.47970 63771 01901 40667 


-0.488 065 354 215 765 843 


-0.49676 28898 97946 30786 


3840 


-0.47970 63771 01901 40689 


-0.488 065 354 215 765 860 


-0.49676 28898 97946 30823 


OO 


-0.47970 63771 01901 4080(3) 


-0.488 065 354 215 765 98(4) 


-0.49676 28898 97946 314(3) 
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TABLE II: The total energies (E) of the bound P{L = 1)— states in the symmetric muonic molec- 
ular ions in muon-atomic units (m^ = l,?i = l,e = l). N designates the number of basis functions 
used in Eq.(2). 





E{ppn) 


E{ddiJ.) 


E{ttiJ.) 


3300 


-0.468 458 430 358 808 027 03 


-0.513 623 952 704 526 3269 


-0.533 263 445 209 533 2938 


3500 


-0.468 458 430 358 808 031 55 


-0.513 623 952 704 526 3342 


-0.533 263 445 209 533 3070 


3700 


-0.468 458 430 358 808 035 14 


-0.513 623 952 704 526 3387 


-0.533 263 445 209 533 3189 


3840 


-0.468 458 430 358 808 037 51 


-0.513 623 952 704 526 3407 


-0.533 263 445 209 533 3254 


oo 


-0.468 458 430 358 808 045(3) 


-0.513 623 952 704 526 39(3) 


-0.533 263 445 209 533 38(3) 


N 


Eipd^^) 


E{pt^i) 


E{dtn) 


3300 


-0.490 664 164 603 504 64 


-0.499 492 024 990 190 10 


-0.523 191 452 003 587 60 


3500 


-0.490 664 164 603 507 57 


-0.499 492 024 990 191 53 


-0.523 191 452 003 588 67 


3700 


-0.490 664 164 603 510 41 


-0.499 492 024 990 192 53 


-0.523 191 452 003 589 54 


3840 


-0.490 664 164 603 511 76 


-0.499 492 024 990 192 96 


-0.523 191 452 003 590 17 


oo 


-0.490 664 164 603 515(1) 


-0.499 492 024 990 195(1) 


-0.523 191 452 003 593(1) 


N 


E*{ddfi) 


E*{dtfi) 


E*{ttn) 


3300 


-0.473 686 731 121 137 629 


-0.481 991 527 054 2451 


-0.489 908 663 057 013 5630 


3500 


-0.473 686 731 121 137 765 


-0.481 991 527 054 3644 


-0.489 908 663 057 013 6844 


3700 


-0.473 686 731 121 137 946 


-0.481 991 527 054 4505 


-0.489 908 663 057 013 8018 


3840 


-0.473 686 731 121 138 063 


-0.481 991 527 054 4894 


-0.489 908 663 057 013 8571 


oo 


-0.473 686 731 121 138 5(3) 


-0.481 991 527 054 9(3) 


-0.489 908 663 057 014 5(2) 
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TABLE III: The total energies (E) of the bound D(L = 2)— states of the the ddiJ,,ttfi and dtij, 
muonic molecular ions in muon atomic units. A'^ designates the number of basis functions used in 
Eq.(2). 





E{ddn) 


E{ttii) 


N 


E{dtii) 


2800 


-0.488 708 327 4382 


-0.512 568 647 4651 


3900 


-0.500 118 078 7334 


3000 


-0.488 708 327 5083 


-0.512 568 647 6583 


4200 


-0.500 118 078 7337 


3200 


-0.488 708 327 5460 


-0.512 568 647 8022 


4500 


-0.500 118 078 7338 


3400 


-0.488 708 327 5584 


-0.512 568 647 9152 


4800 


-0.500 118 078 7338 


3600 


-0.488 708 327 5715 


-0.512 568 648 0085 


5100 


-0.500 118 078 7338 


oo 


-0.488 708 327 598(3) 


-0.512 568 648 210(5) 


oo 


-0.500 118 078 7338(1) 



TABLE IV: The bound state properties X computed for the ground S{L = 0)— state and excited 
S*{L = 0)— state in the pdfj, and dtfj, muonic molecular ions (in muon-atomic units). 



X 


pdji{S{L = 0)-state) 


dttJ,{S{L = 0) -state) 


dtn{S*{L = 0) -state) 




0.6411463600638491(1) 


0.7227000026976390(3) 


0.5146887255965181(3) 




0.7533736114443716(1) 


0.7583156054974041(3) 


0.7053753065541911(3) 


{r2i) 


0.3690963865448171(3) 


0.4038256647760819(5) 


0.2439333237191783(6) 


{rai) 


2.451487643610344(2) 


2.117912271227347(3) 


3.933236044506724(3) 


{rz2) 


2.087699160470998(2) 


2.023720516217468(3) 


2.738751054881910(3) 




3.100710462458351(3) 


2.747914171742117(5) 


5.161229304527515(5) 


(rl) 


8.033494559453071(3) 


5.881854047519130(4) 


22.39719714570313(5) 


(rl) 


10.829021567566262(4) 


8.2873255690653121(6) 


30.631304593819174(6) 




20.6547098820345(2) 


17.4696971262928(3) 


65.2551287149941(3) 


(rli) 


42.147966883621(4) 


27.208344987497(5) 


201.45182645373(6) 


(-^Vf) 0.2806191821887450(7) 


0.3910764626134503(9) 


0.3836327048408983(9) 




0.3674608971218589(7) 


0.4218837688718115(9) 


0.4476331876944995(9) 


(-iVi) 0.46041133009768349(4) 0.50069529322880649(5) 0.50035758500774130(6) 
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TABLE V: The convergence of the {S32), (<^3i) c^nd {S31) expectation values for the ground (bound) 
S{L = 0)— state of the pdfj, molecular ion (in muon-atomic units). 



N (J32) (^31) ^21) 

3300 1.73456203087-10-1 1.17709732798-10-^ 1.46169407-10-^ 
3500 1.73456202965-10-1 1.17709733128-10-1 1.46169370-10"^ 
3700 1.73456202754-10-1 1.17709733014-10-1 1.46169377-10"^ 
3840 1.73456202768-10-1 1.17709733160-10-1 1.46169383-10-^ 



TABLE VI: The levels of hyperfine structure e and hyperfine structure splittings A in the bound 
S{L = 0)— states of the pd/j. and dt/j. ions (in MHz). 



ej=2{pdijL) 1.2519350851-10^ 

ej=i{pdfM) 9.3058194294-10^ A(J = 2 J = 1) 3.2135314217-10^ 

ej=o{pdn) -2.1222395094-10'^ A(J = 1 ^ J = 0) 3.0528214524-10'^ 

ej=\{pdiJL) -2.3097272483-10^ A(J = ^ J = 1) 1. 8748773889- 10^ 

ej=2{dtn) 1.8919590437-10^ 

ej=i{dtii) 1. 5985479092- 10^ A(J = 2 J = 1) 2.9341113453-10^ 

ej^oidtfi) -3.4439378258- lO'^ A(J = 1 ^ J = 0) 5.0424857350-10'^ 

ej=i{dtii) -3.6038337067-10^ A(J = ^ J = 1) 1.5989588085-10^ 



ej=2{dtfi)* 1.8609555434-10^ 

ej=i{dtn)* 1.6554331952-10^ A(J 
ej=o{dtnY -3.4859818025-10^ A(J 
ej=i{dtiiy -3.5950318334-10^ A(J 



= 2 J = 1) 2.0552234821 -10^ 
= 1 J = 0) 5.1414149977-10^ 
= ^ J = 1) 1.0905003094-10^ 
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TABLE VII: The levels of hyperfine structure e and hyperfine structure splittings A in the ground 
S{L = 0)-state of the pt/i ion (in MHz). 



ej_3{ptiJL) 1.5079820356-10'^ 

ej^iiptji) 1.6796717260-106 A(J = | ^ J = 1) 1.3400148630-10^ 
ej^iiptu) -3.1839312439-10^ A(J = i ^ J = i) 3.3518984165-10^ 



TABLE VIII: The convergence of the {dp^) and [5^^) expectation values for the ground (bound) 
S{L = 0)— state of the pp^ molecular ion (in muon-atomic units). 

N ((5,31) (^2l) 

3300 1.315008614364-10^1 3.9370034861-10-5 

3500 1.315008614369-10-1 3.9370034722-10-^ 

3700 1.315008614374-10-1 3.9370034782-10"^ 

3840 1.315008614378-10-1 3.9370034773- 10"^ 



TABLE IX: The hyperfine structure and hyperfine structure splitting of the bound S{L = 0)— states 
of the pp^i and tt^i ions (in MHz) . 







ppii 


ttji 


{ttiiY 


€ T_3 

•^ — 2 




1.256448515-10^ 


1.736310113-10^ 


1.510018118-10^ 






1.772596177-10^ 


1.105911127-101 


1.243070661-101 






-2.513074289-10^ 


-3.472621366-10^ 


-3.020037480-10^ 


A(i- 




1.256271251-10'^ 


1.735309024-10'^ 


1.510016875-10'^ 






2.513251549-10^ 


3.472622472-10'^ 


3.020038723-10^ 
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TABLE X: The hyperfine structure and hyperfine structure splitting of the bound S{L = 0)— states 
of the ddfi ion (in MHz). 



ddjj, (dd^)* 
4.656669271-106 4.023227167-106 
2.328339810-106 2.011617137-10^ 
6.899579465-10° 4.737788941-10° 
-4.656669271-106 -4.023227167-106 
-6.985012530-106 -6.034846672-106 

A(| ^ I) 2.328329461-106 2.01161003-106 

A(| ^ ^) 2.328332911-106 2.01161239-106 

A(i ^ i) 4.656676170-106 4.02323191-106 

A(i -)■ I) 2.328343259-106 2.01161951-106 



e T_6 

J — 2 
•^ — 2 
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